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Introduction 

This manual describes the use of the ANSCARES program to prepare a neutral file of 
FEM stress results taken from ANSYS Release 5.0, 5 in the format needed by the CARES/LIFE 1 
ceramics reliability program. It is intended for use by experienced users of ANSYS and CARES. 
Knowledge of compiling and linking FORTRAN programs is also required. ANSCARES has 
been jointly developed by NASA’s Lewis Research Center and The Advanced Manufacturing 
Center at Cleveland State University. 

Maximum use is made of existing routines (from other CARES interface programs and 
ANSYS routines) to extract the finite element results and prepare the neutral file for input to the 
reliability analysis. FORTRAN and machine language routines as described in Reference 6 are 
used to read the ANSYS results file. Sub-element stresses are computed and written to a neutral 
file using FORTRAN subroutines which are nearly identical to those used in the NASCARES 
(MSC/NASTRAN to CARES) interface. 


* Resident Research Associate - NASA Lewis Research Center 
f Advanced Manufacturing Center 
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Ceramics Reliability Analysis: 

Increasing numbers of components are being made of ceramic materials in order to take 
advantage of higher temperature and corrosion resistance characteristics. Unfortunately, ceramic 
materials do not follow the same simple strength relations that characterize ductile metals. The 
CARES/LIFE (Ceramic Analysis and Reliability Evaluation of Structures Life Prediction 
Program) 1 program has been developed by NASA Lewis Research Center to predict the strength 
and reliability of ceramic components. 

Probabilistic component design involves predicting the probability of failure for a 

f 

thermomechanically loaded component from specimen rupture data. Typically these experiments 
are performed using many simple geometry flexural or tensile test specimens. A static, dynamic, 
or cyclic load is applied to each specimen until fracture. Statistical strength and SCG (fatigue) 
parameters are then determined from these data. Using these parameters and the results obtained 
from a finite element analysis, the time-dependent reliability for a complex component geometry 
and loading is then predicted. Appropriate design changes are made until an acceptable 
probability of failure has been reached. This design methodology combines the statistical nature 
of strength-controlling flaws with the mechanics of crack growth to allow for multiaxial stress 
states, concurrent (simultaneously occurring) flaw populations, and subcritical crack growth. 
These issues are addressed within the CARES/LIFE program. 

The stress distribution within the component is provided by a file extracted from the 
results of a finite element method (FEM) calculation. Interfaces have been developed to the 
co mm ercial FEM programs ABAQUS, Version 4.9-5.X (ABACARES) 2 , MSC/NASTRAN 
version 67 (NASCARES) 3 , and ANSYS, Release 4.4 (ANSYS/CARES) 4 . The ABAQUS and 
NASTRAN versions prepare a neutral file, in standardized format, which is accepted by the most 
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recent reformulated CARES/LIFE program, while the ANSYS version uses the older concept of 
reading the FEM results through subroutines called from inside the CARES program. This 
current effort develops an interface that is compatible with the most recent versions of ANSYS 
and CARES/LIFE. A translator for ANSYS version 4.4A to CARES/LIFE is also available. 

Neutral Data Base Concept 

The concept of a formatted neutral data base allows for the interfacing of CARES/LIFE 
to several finite element packages. The results from the finite element analysis (available in the 
form of standard output or binary files and/or plot data files) have to be interpreted by an 
interface program. Due to the different formats used for output of the results in each finite 
element program, this interpreter program has to be adapted to the finite element software used. 

The finite element results required for the reliability evaluation are assembled into the 
formatted neutral data base. A further advantage of this approach is the ease of transfer of this 
formatted neutral database (ASCII file) to different computer systems. Interpreter control 
parameters, such as the number of element groups to be considered in the reliability analysis, 
flags for accounting for temperature field gradients, etc., are usually input via an interpreter 
control file. 

A subelement technique has been implemented. The reliability analysis is now performed 
at each Gauss integration point instead of using element averaging techniques. In the context of 
finite element analysis, stresses are determined at the Gaussian integration points where the local 
stiffness matrix has been evaluated. The subelement technique implies that each Gaussian 
integration point corresponds to a subelement. The subelement volume or area is defined as the 
contribution of the integration point to the element volume in the course of the numerical 
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integration procedure. 

The location of the Gaussian integration point in the natural space of the finite element, 
as well as corresponding weight functions have to be considered when the subelement volume 
is calculated. For this reason, the number of subelements in each element depends on the 
integration order chosen, and as a consequence, the element type. The increased number of 
points for stress evaluation accounts for the variation in stress over the element. Thus, 
considerable improvements in the accuracy of the reliability analysis have been realized. 3 

The neutral data base supplies the finite element results and data required for the 
reliability analysis, including the stress tensors and temperatures at the integration points. The 
structure of the neutral data base is optimized with respect to memory. The finite element data 
is arranged within the neutral data base using the following hierarchy: element groups, elements 
and subelements. The element group data contain information regarding the number of elements 
wi thin the group. In addition to the element types available in the previous version of CARES 
(three-dimensional brick elements for volume flaw analysis, membrane elements for surface flaw 
analysis), two-dimensional finite elements are also implemented. By assuming plane stress 
conditions or exploiting symmetry of the structure, the size of the finite element model can be 
reduced using a two-dimensional mesh. Thus, the number of degrees of freedom could be 
substantially reduced. 

The two-dimensional modeling options available (in plane loading vs. bending) depend 
on the method used to determine the volume. For ANSYS 5.0, the volume is calculated by 
multiplying the area of a subelement by the thickness of the element. Therefore, bending effects 
are not taken into account because the stress is assumed to be constant through the thickness. 

Information pertaining to the elements include the number of subelements and the material 
identification n um ber. Finally, the subelement group data contain information regarding 
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subvolume, stress tensor (i.e., the stress tensor at the integration point), and subelement 
temperature. The stresses and temperatures are assumed to be constant over the subelement 
considered. 

Based on the global hierarchy outlined above, the neutral file data is arranged into records 
with standard FORTRAN formats: the first record (format A80) specifies the title for the 
particular reliability evaluation. The next entry (format 215) contains the number of element 
groups (both volume as well as surface (shell) elements) to be considered in the reliability 
analysis and the a flag indicating combined volume as well as surface flaw reliability analysis for 
shell elements. 

For every element groups, a loop is started with the next statement (format 415): the 
identification number of the element group considered, the element type corresponding to 
PATRAN element type codes, the number of elements within this element group, and a shell 
element type indicator have to be input. If shell elements are used in the analysis, the indicator 
has to be set to 1 . Information pertaining to the element level is read in a loop over the number 
of elements within the element group considered (format 4I5,3F15.0). 

The element level data contain the identification number of the element considered, the 
number of subelements within this element, the material identification number for volume flaw 
analysis, the material identification number for surface flaw analysis, the element volume, the 
element thickness only in case of shell elements, and the averaged temperature of this element. 
The identification number of the element should be the same as the corresponding element 
number in the finite element mesh, if further postprocessing of the reliability results are planned. 
The number of subelements within this element depends on the Gaussian integration order chosen 
in the finite element analysis and on the element type. If stress data is only available at the 
center of the element (standard case in NASTRAN analyses prepared for previous CARES 
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versions), one subelement is specified. The material identification numbers must correspond to 
the material numbers used in the finite element analysis as well as the material numbers on the 
control file (containing the Weibull parameters). Element volume and element temperature are 
necessary only for element group information summaries in the printout file. For the reliability 
evaluation the subelement data are directly used. 

A third loop is started for reading the subelement level data. For volume type elements 
the identification number of the subelement, the subelement volume, the subelement temperature 
and the full stress tensor at the integration point (defined by a„, a^, and in the 

local coordinate system) have to be specified (format I5,5F15.0,/,35X,3F15.0). 

For shell type elements plane stress conditions are assumed. For that reason only o„, 
and Qxj, tire necessary for the definition of the subelement stress tensor. Instead of the subelement 
volume, the subelement area has to be specified (I5,5F15.0). 

Typically, subelement volumes are not included with standard finite element output. Thus, 
the volume of each subelement (corresponding to a Gauss integration point) is calculated in the 
interpreter program using the shape functions inherent to the element type. In the usual context 
of finite element methods, the volume of a three-dimensional element, (i. e., brick, wedge, 
pyramid or tetrahedron) is calculated after transformation into the natural coordinate space, 

1 1 1 

V= f f fdetJ (rj,t) dr ds dt (1) 

-l -l -l 

where J is the Jacobian operator and r,s,t are the natural coordinates. By applying a Gaussian 
integration scheme, the integral can be expressed as 
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( 2 ) 


v = E E E d “ J(w>> w, ^ >n 

i-l j-1 k-1 

where n is the integration order, r;,Sj,t k are the coordinates of the Gaussian integration points in 
natural space and the W’s are the weight functions. The area of an element (i. e., quadrilateral 
or triangle) is determined in a similar manner 

D D 

A= £ £ det Kr.jp W { Wj (3) 

i-i j-i 

The volume of an axisymmetric element can also be determined similarly 

V= R £ £ det J (V0 W t Wj (4) 

i-i j-i 

where R is the radial coordinate of the integration point in the global coordinate system. 

ANSYS/CARES Analysis Procedure: 

Fracture probability prediction on a ceramic part begins with the execution of the ANSYS 
program using a finite element model of the part in question. This model will typically consist 
of volume elements with shell elements attached at the surfaces to provide both volume and 
surface data to CARES/LIFE. Supported element types are listed in Table 1. Shell elements 
should have membrane properties only and thicknesses near zero to avoid affecting the stiffness 
of the component. It is vital that the material identifications used in the ANSYS element 
descriptions agree exactly with those contained in the CARES/LIFE template file (MATID 
keyword in the Material Control Input). 

Standard results utilized from the ANSYS linear elastic analysis results file include nodal 
coordinates, temperatures, and stress components at all nodes. The ANSCARES program is then 
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used to compute subelement stresses, volumes, and temperatures as needed by the CARES/LIFE 
program. All needed data are read from the ANSYS results file (FILE.RST), and no user input 
is required except to change file naming and analysis type defaults, if desired. ANSCARES 
interprets one load case per run, therefore ANSCARES does not handle multiple load cases. 

The final step is to execute the CARES/LIFE program using the neutral file prepared by 
ANSCARES and a template file containing material fracture properties. The template file also 
controls the type of analysis and the results computed. Figure 1 shows a flow diagram of a 
typical reliability analysis using finite element stress and volume results. 

« 

Installing ANSCARES: 

The ANSCARES program is provided on disk with an example problem. ANSCARES 
consists of a main program and nine subroutines, all supplied in FORTRAN 77 source code in 
two modules called ANSCARES 1. FOR and ANSCARES2.FOR (anscaresl.f and anscares2.f: 
renamed for unix systems). This module should be copied into a suitable directory and compiled. 

ANSCARES makes use of system-dependent, machine language subroutines included with 
the ANSYS distribution tape. These routines are designed to read data from the ANSYS results 
file. The specific subroutines needed are BINCLO, BININI, BINIQR, BINRD, BINSET, 
EXINC4, and SYSIQR, all of which are included in the machine language file BINLIB.OLB. 
This file must be copied from the ANSYS distribution tape to the directory used to run the 
ANSCARES program. 

These machine language ANSYS routines must be linked with the compiled FORTRAN 
object code to create an executable version of ANSCARES. If they are all in the same library, 
the command to do this on a VAX/VMS computer is 
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LINK ANSCARES,BINLIB/LIBRARY 


for a silicon graphics and possibly other unix systems, the command is 

f77 -mp -g -o anscares anscaresl.o anscares2.o /ansys50/objects/binlib.a 

The options for the f77 command are those which were used to link the program on a SGI 
Personal IRIS. Linking operations are computer specific, and can be expected to make use of 
a different command on other computer systems. However, these are described in the appropriate 
ANSYS INSTALLATION GUIDE for your system under Stand-Alone Utility Programs, BINLIB 
Subroutines. 

The ANSCARES disk also contains the ANSYS input file for the example problem, 
below, along with the corresponding CARES/LIFE template file. To confirm the proper 
installation, ANSYS, ANSCARES, and CARES/LIFE should be executed in turn using these 
files. Copies of the resulting ANSCARES output neutral file (ANSCARES.NEU) and 
CARES/LIFE output (CARES. OUT) are included for comparison. 

ANSYS considerations: 

At present, analysis employing eight ANSYS element types has been implemented into 
the ANSCARES program. These element types allow three-dimensional and axisymmetric 
analyses. Table 1 describes the ANSYS Release 5.0 element types supported by the ANSCARES 
program. To provide the stress results needed by ANSCARES, all element types must have the 
element description variable KEYOPT(5) set to a value of two (2). Element types that use the 
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variable KEY0PT(6) should have its value remain at the default value of zero (0). 

ANSCARES program execution 

The ANSCARES program will interactively prompt for the names of the ANSYS stress 
and geometry results file to be read and the CARES neutral file to be generated. Depressing the 
ENTER key for each prompt will accept the default names of "FILE.RST" and "CARES.NEU". 
Additional prompts are for element or subelement analysis and waming/progress messages to be 

sent to the monitor or to a file. Subelement analysis is highly recommended. Again pressing the 

« 

ENTER key will accept the default values of subelement analysis (ICARES=2) and display of 
messages on the monitor. A response of "n" to the prompt will allow the user to specify that 
messages are saved in the "ANSCARES. OUT" file instead, when IOUT is set to 9 instead of 
(default) 6. Depending on the computer system, it may be necessary to change the default value 
for monitor display to some value different from 6. 

There is no specific limit on the number of nodes and elements permitted in the finite 
element model. All needed node and element data is stored in a general purpose array (the A 
array). The maximum number of nodes and elements is limited only by the amount of space 
available in this array. Of course, some element types require more data for each element and 
use more space than others. The declared size of the array can be increased or decreased as 
necessary. The maximum size of the A array is set by the parameter MTOT in the ANSCARES 
program (main routine) and is related to the amount of computer resources available. After 
changing the MTOT parameter, the ANSCARES main program, but not necessarily the 
subroutines, must be re-compiled and all programs re-linked to establish the new limits. 

The parameter ITOT in the ANSCARES program (main program) is used to define the 
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size of the temporary buffer used when reading the ANSYS results file. The ANSYS results 
variable RECLNG contains the minimum size of buffer needed. ITOT is currently set to the size 
needed by the VAX computer version of ANSYS (7680 words). The ANSCARES program will 
print a warning message and exit if the computer system in use requires more space. To 
minimiy^ computer storage it is necessary to determine the minimum buffer needed. This can 
be accomplished by setting the ITOT variable to 1 which will result in a message giving the 
required size when ANSCARES is executed. ITOT can then be given that value. As with 
MTOT, after changing the ITOT parameter the ANSCARES main program must be re-compiled 
and all programs re-linked to establish the new buffer size. 

ANSCARES Program Description: 

Machine language routines supplied with the ANSYS program are called as needed to read 
the ANSYS analysis results file. This file is described in Chapter 6 of Ref [6]. The routines 
used to access this file are documented in Chapter 5, also in Ref [6]. 

The ANSCARES main program serves to establish sizes for arrays, call the subroutines 
that read the ANSYS results file, and call the subroutines that prepare the CARES/LIFE neutral 
file. Figure 2 shows the major functions of the various modules in the ANSCARES program. 
ANSYS results are read by the subroutines EANSYS, RDGEOM, and RDSOLU. These three 
routines are based on and remain very similar to the example programs beginning on page A-21 
of Ref. [6]. 

The EANSYS subroutine opens the ANSYS results file (FILE.RST) and determines the 
number of nodes and elements in the model, and other parameters. Subroutine RDGEOM is then 
called by EANSYS to extract coordinate locations, element connectivity, and properties. This 
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data is stored in the A array for later use. A count is also made of the various kinds of elements 
in the model. EANSYS then calls subroutine RDSOLU to extract stress and temperature results, 
which are also stored in the A array. This completes reading of the data from the ANSYS results 
file. 

The HEXA, PENTA, TETRA, QUAD8, TRIA6, and TRIAX subroutines are then called 
by ANSCARES as needed, depending on the kinds of elements in the model. These subroutines 
are very similar to the same-named subroutines in the previous NASCARES program. 3 

« 

Example Problem: 

Sample problems were executed using the ANSYS-ANSCARES-CARES/LIFE sequence 
to verify the accuracy of the analysis as compared to interfaces for earlier versions of ANSYS. 
Models including elements with and without mid-side nodes and using solid, shell, and 
axisymmetric element types confirmed program operation. One of the models, EXAMPLE 1 from 
Reference 4, is described below. 

EXAMPLE 1 is a transversely loaded ceramic circular plate. In ANSYS this is treated as 
a sector suitably restrained to preserve symmetry. The plate is subdivided into 48 equal sectors, 
making each sector 7.5° wide. The model of the sector is divided into 40 three-dimensional solid 
elements and 10 shell elements, all with mid-side nodes. The four solid elements and the shell 
element at the center of the model use duplicated node numbers to make five sided prism solid 
and triangular shell elements. The remaining 36 solid and 9 shell elements remain as solids with 
six sides and shells with four edges. 

The circular plate has an outside diameter of 2.00 inches and a thickness of 0.071 inches. 
Figure 3 is a side view of the sector, showing the element numbers. The four layers shown are 
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made up of the solid elements, with the ten shell elements attached to the bottom of the lowest 
layer. Figures 4 and 5 show additional details of the node and element numbering. The sector 
is simply supported between elements 9 and 10 (equivalent to a ring support at a radius of 0.92 
inches). A uniform pressure of 220 psi is applied to elements 31 through 39. Temperature is 
a uniform 70°F in all elements. Material properties used in this problem are given below: 

- Material Number = 300 (Must agree with CARES template) 

- Young’s Modulus = 5.87 x 10 7 psi 

- Density =0.5 pounds/cubic inch 

- Poisson’s Ratio = 0.25 

The entire ANSYS input file is included in Appendix A. 

This example is provided to demonstrate model generation and neutral file design. The 
CARES/LIFE template file is supplied to run CARES/LIFE using the generated neutral file. The 
background for this example is described in reference [7]: Example 1 - Transversely Loaded 
Circular Disk. 
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ANSYS 

Element 

Number 



Type 

of Nodes 

Description 

Comment 

42-PLANE 

4 

Axisymmetric - 4 Edges 

Not supported 


3 

Axisymmetric - 3 Edges 

KEYOPT(3) = 1 

45-SOLID 

8 

Brick - 6 Sides 



6 

Prism - 5 Sides 



4 

Tetrahedron - 4 Sides 


63-SHELL 

4 

Shell - 4 Edges 



3 

Shell - 3 Edges 

KEYOPT(l) = 0 

72-SOLID 

4 

Tetrahedron - 4 Sides 


82-PLANE 

8 

Axisymmetric - 4 Edges 

Not Supported 


6 

Axisymmetric - 3 Edges 

KEYOPT(3) - 1 

92-SOLID 

10 

Tetrahedron - 4 Sides 


93-SHELL 

4 

Shell - 4 Edges 



3 

Shell - 3 Edges 


95-SOLID 

20 

Brick - 6 Sides 



15 

Prism - 5 Sides 



13 

Pyramid - 5 Sides 

Not Supported 


10 

Tetrahedron - 4 Sides 



Note: For all of the ANSYS element types above: 

KEYOPT(5) + 2 
KEYOPT(6) = 0 


Table 1 - ANSYS element types supported by ANSCARES 
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Figure 1 - Block diagram for the analysis and time-dependent reliability evaluation of 
ceramic components. This analysis includes the CARES/LIFE integrated design 
program, finite element analysis and, PATRAN PLUS pre- and post-processing. 
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Begin AN SC ARES program 
Check for ANSYS results file 


Call E ANSYS Subroutine 

Extract disk pointers and model parameters 

Call RDGEOM Subroutine 

Extract and save node coordinates 
For each element: 

Extract and save element properties 
Extract element connectivity 

Save in temporary space 
Add to element type counter 
For each element type: 

Eliminate duplicate element nodes 
Move connectivity to save array 

Call RDSOLU Subroutine 
For each element: 

Extract and save element stresses 
For each element: 

Extract temperatures at nodes 
Save temperatures in node order 

Write neutral file header 

For each element type in the model: 

Call appropriate element type subroutine 

(HEXA, PENT A, TETRA, QUAD8, TRIA6, or TRIAX) 
For each element of this type: 

Compute sub-element stresses 
Compute sub-element volumes 
Compute average temperature 
Write results to neutral file 


Close files 

End AN SC ARES program 


Figure 2 - ANSCARES Program Logic 


Version 5.0 


- 16 - 


May, 1995 



Distance of Solid Elements 
from Support Surface 


.071 


Uniform Pressure of 220 psi 



Simple Support 

Radial Distance from Center of Disc (in.) 


Figure 3 - Side view of a sector of the ceramic plate 


Preliminary Version 


- 17 - 


April, 1994 




















Figure 5 - Element plot of bottom shell layer 
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Appendix A 

ANSYS Input File for Examplel 

/BATCH,LIST 

/PREP7 

/TITLE, STRESS ANALYSIS OF A SIMPLY SUPPORTED DISC, UNIQUE ELEMENT PRESSURE 

ET,1,95 

KEYOPT, 1,5,2 

ET,2,93 

KEYOPT,2,5,2 

EX,300,5.87E07 

DENS,300,0.5 

NUXY,300,0.25 

TREF,70 

csys,l 

N,1,0.,0.,0. 

N,2,0.01,0.,0. 

N,3,0.02,0.,0. 

N,4,0.04,0.,0. 

N,5,0.06,0.,0. 

N,6,0.10,0.,0. 

N,7,0.14,0.,0. 

N,8,0.22,0.,0. 

N,9,0.30,0.,0. 

N,10,0.42,0.,0. 

N,ll,0.54,0.,0. 

N,12,0.66,0.,0. 

N,13,0.78,0.,0. 

N, 14,0.82, 0.,0. 

N,15,0.86,0.,0. 

N,1 6,0.88,0., 0. 

N,17,0.90,0.,0. 

N,1 8,0.9 1,0.,0. 

N,19,0.92,0.,0. 

N,20,0.96,0.,0. 

N.21,1.00,0.,0. 

NGEN.2,72, 1,2 1 , 1 ,,,0.002 
NGEN,2,72,73,93,1,„0.008 
NGEN,2,72, 145, 165,1, ,,0.0255 
NGEN.2,72,2 1 7,237, 1 ,,,0.0355 
N,52,0.,0., 0.001 
N,53,0.02,0., 0.001 
N,54,0.06,0.,0.00 1 
N,55,0.14,0.,0.001 
N,56, 0.30,0., 0.001 
N,57,0.54,0.,0.00 1 
N,58,0.78,0.,0.00 1 
N,59,0.86,0.,0.00 1 
N,60,0.90,0.,0.00 1 
N,61,0.92,0.,0.001 
N,62,1.00,0.,0.001 
NGEN,2,72,52,62, 1 ,,,0.005 
NGEN.2,72, 1 24, 1 34, 1 ,„0.0 1 675 
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NGEN,2,72, 1 96,206, 1 ,,,0.03025 
N,22,0.02, 3.75,0. 

N^3, 0.06, 3.75,0. 

N,24,0.14,3.75,0. 

N,25,0.30,3.75,0. 

N,26,0.54,3.75,0. 

N.27,0.78, 3.75,0. 

N,28,0.86,3.75,0. 

N,29,0.90,3.75,0. 

N,30,0.92,3.75,0. 

N^l, 1.00,3.75,0. 

NGEN.2,72, 22,3 1,1 ,,,0.002 
NGEN,2,72,94, 1 03, 1 ,,,0.008 
NGEN,2,72, 1 66, 1 75, 1,„0.0255 
NGEN ,2,72,238,247,1, ,,0.0355 
NGEN,2,30,2,21,1„7.5, 
NGEN,2,72,32,5 1 , 1 ,,,0.002 
NGEN,2,72, 1 04, 1 23, 1 ,,,0.008 
NGEN,2,72, 1 76, 1 95, 1 ,,,0.0255 
NGEN,2,72,248,267, 1 ,,,0.0355 
NGEN,2,10, 53,62, 1„7.5, 
NGEN,2,72,63,72, 1 ,,,0.005 
NGEN,2,72, 1 3 5, 1 44, 1 ,„0.0 1 675 
NGEN.2, 72,207,2 1 6, 1 ,,,0.03025 
nrot,l,339,l 
TYPE,1 
MAT,300 

E,U, 33,33,73,75, 105, 105 
EMORE.2,22,3 3 ,32,74,94, 1 05, 1 04 
EMORE, 52,53, 63,63 
E,3,5,35,33,75,77, 1 07, 1 05 
EMORE,4,23 ,34,22,76,95, 1 06,94 
EMORE.53,54,64,63 
E,5,7,37,35,77,79, 1 09, 1 07 
EMORE,6,24,36,23,78,96, 1 08,95 
EMORE,54,5 5 ,65 ,64 
E,7,9,39,37,79,81,l 11,109 
EMORE.8,25, 38,24,80, 97,1 10,96 
EMORE,55,56,66,65 
E,9,l 1,41,39,81,83,1 13,1 11 
EMORE, 10,26,40,25,82,98,1 12,97 
EMORE,56, 57, 67,66 
E,1 1,13,43,41,83,85,1 15,1 13 
EMORE, 1 2,27,42,26,84,99, 1 14,98 
EMORE,57,58,68,67 
E, 13, 15,45,43,85, 87,1 17,1 15 
EMORE, 14,28,44,27, 86, 100,1 16,99 
EMORE.58, 59,69, 68 
E,15,17,47,45,87,89,l 19,1 17 
EMORE, 16,29,46,28, 88, 101, 118,100 
EMORE.59, 60,70, 69 
E, 17, 19,49,47,89, 9 1,1 2 1,1 19 
EMORE, 1 8,30,48,29,90, 1 02, 1 20, 1 0 1 
EMORE,60,6 1 ,7 1 ,70 
E, 19,2 1,5 1,49, 91, 93, 123,121 
EMORE,20,3 1 ,50,3 0,92, 1 03 , 1 22, 1 02 
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EMORE,6 1,62,72, 71 

E,73,75,105,105,145,147,177,177 

EMORE,74,94, 1 05, 1 04, 1 46, 1 66, 1 77, 1 76 

EMORE.124, 125, 135,135 

E,75,77,107, 105, 147, 149,179,177 

EMORE,76,95, 1 06,94, 1 48, 1 67, 1 78, 1 66 

EMORE,125,126,136,135 

E,77,79,l 09, 1 07, 1 49, 1 5 1 , 1 8 1 , 1 79 

EMORE,78,96, 1 08,95, 1 50, 1 68, 1 80, 1 67 

EMORE,126,127,137,136 

E,79,81, 11 1,109, 151, 153, 183, 181 

EMORE,80,97, 1 1 0,96, 1 52, 1 69, 1 82, 1 68 

EMORE,127,128,138,137 

E,8 1,83,1 13,1 1 1,153,155,185,183 

EMORE.82,98, 1 1 2,97, 1 54, 1 70, 1 84, 1 69 

EMORE.128,129,139,138 

E,83,85, 1 1 5, 1 1 3, 1 55, 1 57, 1 87, 1 85 

EMORE,84,99, 1 1 4,98, 1 56, 1 7 1 , 1 86, 1 70 

EMORE.129,130,140,139 

E,85,87, 1 1 7,1 1 5, 1 57, 1 59, 1 89, 1 87 

EMORE.86, 1 00, 1 16,99, 1 58, 1 72, 1 88, 1 7 1 

EMORE, 130, 13 1,14 1,140 

E,87,89,l 19,1 17,159,161,191,189 

EMORE,88, 101, 118,100,160,173,190,172 

EMORE,131,132,142,141 

E,89,91, 121,1 19,161,163,193,191 

EMORE,90, 1 02, 1 20, 1 0 1 , 1 62, 1 74, 1 92, 1 73 

EMORE.132,133,143,142 

E,91,93, 123,121, 163, 165,195, 193 

EMORE.92,103, 122, 102, 164, 175, 194,174 

EMORE, 133, 134, 144, 143 

E, 145, 1 47, 1 77, 1 77,2 1 7,2 1 9,249,249 

EMORE, 1 46, 1 66, 1 77, 1 76,2 1 8,23 8,249,248 

EMORE, 1 96, 1 97,207,207 

E, 1 47, 1 49, 1 79, 1 77,2 1 9,22 1 ,25 1 ,249 

EMORE, 148, 167,178,1 66,220,239,250,23 8 

EMORE, 1 97, 1 98,208,207 

E, 149,151, 181, 179, 221,223,253, 251 

EMORE, 1 50, 1 68, 1 80, 1 67,222,240,252,239 

EMORE, 198, 199,209,208 

E, 15 1,1 53, 183, 18 1,223, 225, 255,253 

EMORE, 1 52, 1 69, 1 82, 1 68,224,24 1 ,254,240 

EMORE, 1 99,200,2 1 0,209 

E, 153, 155, 185, 183, 225, 227, 257,255 

EMORE, 1 54, 1 70, 1 84, 1 69,226,242,256,24 1 

EMORE,200, 20 1,2 11,210 

E,155,157, 187, 185,227,229,259,257 

EMORE, 1 56, 1 7 1 , 1 86, 1 70,228,243,258,242 

EMORE.20 1,202,2 12,2 11 

E, 1 57, 1 59, 1 89, 1 87,229,23 1 ,26 1 ,259 

EMORE,l 58, 172,1 88,1 7 1,230,244,260,243 

EMORE,202,203 ,2 1 3,2 1 2 

E, 159, 161, 19 1,1 89,23 1,233,263,261 

EMORE, 1 60,1 73, 1 90, 1 72,232,245,262,244 

EMORE,203,204,2 1 4,2 1 3 

E, 16 1,163, 193, 19 1,233,23 5, 265,263 
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EMORE, 1 62, 1 74, 1 92, 1 73 ,234,246,264,245 

EMORE.204,205 ,215,214 

E, 163,165, 195, 193,235,237,267,265 

EMORE, 1 64, 1 75 , 1 94, 1 74,236,247,266,246 

EMORE,205,206,2 16^215 

E2 1 72 1 9^49,249,289,29 1 ,32 1 ,32 1 

EMORE2 1 8,23 8,249248290,3 1 0,32 1 ,320 

EMORE268269,279279 

E2 1 922 1 25 1 ,249,29 1 293 ,323,32 1 

EMORE220239,250238292,3 1 1,322,3 10 

EMORE269270.280279 

E22 1223253 ,25 1293,295,325,323 

EMORE222,240,252239294,3 12,324,31 1 

EMORE,270,27 1,28 1,280 

E223225,255,253295,297,327,325 

EMORE224.24 1 ,254,240296,3 13,326,3 12 

EMORE,271272,282281 

E,225227,257,255297,299,329,327 

EMORE,226242,256,24 1 ,298,3 1 4,328,3 1 3 

EMORE.272273,283282 

E227229,259257,299,30 1 ,3 3 1 ,329 

EMORE,228243,258,242,300, 3 15,330,3 14 

EMORE.273,274284283 

E229.23 1,26 1,259,30 1,303, 333,331 

EMORE23 0,244,260243,302,3 1 6,3 32,3 1 5 

EMORE274,275,285284 

E231233263.261, 303, 305, 335,333 

EMORE232245,262,244,304,3 17,334,3 16 

EMORE275276,286,285 

E233235265,263,305,307,337,335 

EMORE234,246,264,245, 306, 3 18,336,3 17 

EMORE276,277,287,286 

E,235237,267, 265, 307, 309,339, 337 

EMORE,236,247,266,246, 308,3 1 9,338,3 1 8 

EMORE,277,278,288287 

TYPE2 
MAT, 300 
R, 1,0.000001 
E,l, 3,33,33,222, 33, 32 
E,3, 5,35,33,423, 34,22 
E,5, 7,37,35,624, 36,23 
E,7, 9,39,37,8,25, 3824 
E,9, 11,41,39,10,26,4025 
E,1 1,13,43,41,12,27,4226 
E, 13,15,45,43,14,28, 44,27 
E,15, 17,47,45, 16,29, 46,28 
E,17,19,49,47,l 8,30,4829 
E,1921,51,49,20,3 1,50,30 
WSTART,1,12 
WMORE2220 
WMORE,32,32 
WMORE,32,32,l 
WSTART.3 ,322,30 
WMORE, 4,4,30 
WMORE22,22,l 
WSTART, 5,5,2,2,30 
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WMOR£,6,6,30 
WMORE.23,23,1 
WSTART.7, 7,2,2,30 
WMORE.8,8,30 
WMORE,24,24,l 
WSTART.9,9,2,2,30 
WMORE,l 0,1 0,30 
WMORE.25,25,1 
WSTART, 11,11 ,2,2,30 
WMORE, 12, 12,30 
WMORE.26,26,1 
WSTART, 13, 13, 2 ,2 ,30 
WMORE, 14, 14, 30 
WMORE,27,27, 1 
WSTART, 15, 15, 2,2, 30 
WMORE, 16, 16,30 
WMORE,28,28,l 
WSTART, 17, 17,2,2,30 
WMORE, 18, 18,30 
WMORE.29,29,1 
WSTART, 19, 19,2,2,30 
WMORE, 20,20,30 
WMORE, 30, 30,1 
WAVES 
NLIST.ALL 
ELIST.ALL 

/VffiW„-l,2,-l $ENUM,1 SEPLOT 

cp,l,uz^,32 

cp,2,uz,3,22,33 

CP,4,UZ,4,34 

CP,5,UZ,5,23,35 

CP,7,UZ,6,36 

CP,8,UZ,7,24,37 

CP,10,UZ,8,38 

CP,1 1,UZ,9,25,39 

CP, 13,UZ, 10,40 

CP, 14,UZ,1 1,26,41 

CP, 16,UZ, 12,42 

CP, 1 7,UZ, 13^27,43 

CP, 19,UZ, 14,44 

CP,20,UZ,1 5,28,45 

CP,22,UZ, 16,46 

CP,23,UZ,17,29,47 

CP.25.UZ, 18,48 

CP,26,UZ,20,50 

CP,27,UZ,2],31,51 

cp,29,uX,2,32 

cp,30,uX,3 ,22,33 

CP,32,UX,4,34 

CP,33,UX,5,23,35 

CP,35,UX,6,36 

CP,36,UX,7,24,37 

CP,38,UX,8,38 

CP,39,UX,9,25,39 

CP,41,UX,10,40 

CP,42,UX,1 1,26,41 
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CP.44.UX, 12,42 
CP.45.UX, 1 3 ,27,43 
CP,47,UX, 14,44 
CP,48,UX,1 5,28,45 
CP,50,UX, 16,46 
CP,51,UX,17,29,47 
CP,53,UX, 18,48 
CP,54,UX,20,50 
CP,55,UX,2 1,3 1,51 
CP,56,UZ,74,104 
CP,57,UZ,75,94,105 
CP,58,UZ,76,106 
CP,59,UZ,77,95,107 
CP,60,UZ,78,108 
CP, 6 1 ,UZ,79,96, 1 09 
CP,62,UZ,80,1 10 
CP,63,UZ,8 1,97,1 11 
CP,64,UZ,82,1 12 
CP,65,UZ,83,98,1 13 
CP,66,UZ,84,1 14 
CP,67,UZ,85,99,1 15 
CP,68,UZ,86,1 16 
CP,69,UZ,87, 1 00, 1 1 7 
CP,70,UZ,88, 1 18 
CP,7 1,UZ,89,101,119 
CP,72,UZ,90,120 
CP,73,UZ,91, 102,121 
CP,74,UZ,92, 122 
CP,75,UZ,93, 1 03, 1 23 
CP,76,UX,74, 1 04 
CP,77,UX,75,94, 1 05 
CP,78,UX,76,106 
CP,79,UX,77,95, 1 07 
CP,80,UX,78,108 
CP,8 1 ,UX,79,96, 1 09 
CP,82,UX,80,1 10 
CP, 83,UX,8 1,97,1 11 
CP,84,UX,82,1 12 
CP,85,UX,83,98,1 13 
CP,86,UX,84,1 14 
CP,87,UX,85,99,1 15 
CP,88,UX,86,1 16 
CP,89,UX, 87,1 00,1 17 
CP,90,UX,88,1 18 
CP,91,UX,89,101,1 19 
CP,92,UX,90,120 
CP,93,UX,91,102,121 
CP,94,UX,92, 1 22 
CP,95,UX,93 , 1 03 , 1 23 
CP,96,UZ, 1 46, 1 76 
CP,97,UZ, 1 47, 1 66, 1 77 
CP,98,UZ,148,178 
CP,99,UZ, 1 49, 1 67, 1 79 
CP, 1 00,UZ, 150,180 
CP,101,UZ,15 1,168, 181 
CP,1 02,UZ, 1 52, 1 82 
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CP, 103,UZ, 153,1 69, 1 83 
CP, 104,UZ, 154,184 
CP,105,UZ,155,170,185 
CP,106,UZ,156,186 
CP, 1 07,UZ, 1 57, 1 7 1 , 1 87 
CP,108,UZ,158,188 
CP, 1 09, UZ, 1 59, 1 72, 1 89 
CP, 1 1 0,UZ, 1 60, 1 90 
CP,1 1 1,UZ,161,173,191 
CP, 1 1 2,UZ, 1 62, 1 92 
CP,1 13,UZ,163,174,193 
CP,114,UZ,164,194 
CP,1 15, UZ, 165, 175, 195 
CP,1 16, UX, 146, 176 
CP,1 17,UX,147,166,177 
CP,U8,UX, 148,178 
CP, 1 1 9,UX, 1 49, 1 67, 1 79 
CP, 120,UX, 1 50, 1 80 
CP,121,UX,1 5 1,168,1 81 
CP, 1 22, UX, 1 52, 1 82 
CP,123,UX,153,169,183 
CP,124,UX, 154,184 
CP,125,UX,155,170,185 
CP,1 26,UX, 1 56, 1 86 
CP,127,UX, 157,171, 187 
CP,128,UX,158,188 
CP,1 29,UX, 1 59, 1 72, 1 89 
CP,130,UX, 160,1 90 
CP,131,UX,161,173,191 
CP,1 32,UX, 1 62, 1 92 
CP,133,UX,163,174,193 
CP,134,UX, 164,1 94 
CP,135,UX,165,175,195 
CP,136,UZ,2 18,248 
CP,1 37,UZ,2 1 9,23 8,249 
CP,138,UZ,220,250 
CP,139,UZ, 22 1,239,251 
CP, 1 40,UZ,222,2 5 2 
CP, 141 ,UZ,223 ,240,253 
CP, 1 42,UZ,224,254 
CP, 1 43 ,UZ,225,24 1 ,25 5 
CP,144,UZ,226,256 
CP,145,UZ,227,242,257 
CP,146,UZ,228,258 
CP,147,UZ, 229,243, 259 
CP,148,UZ,230,260 
CP, 149,UZ,23 1 ,244,26 1 
CP, 1 50,UZ,23 2,262 
CP, 1 5 1 ,UZ,233,245 ,263 
CP, 1 52,UZ,234,264 
CP, 153 ,UZ,23 5,246,265 
CP,154,UZ,236,266 
CP,155,UZ,237,247,267 
CP, 1 56,UX,2 1 8 ,24 8 
CP,157,UX,219, 238,249 
CP,158,UX,220,250 
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CP, 1 59,UX,22 1 ,239,25 1 

CP, 1 60,UX,222,252 

CP,161,UX,223,240,253 

CP,162,UX,224,254 

CP,163,UX,225,241,255 

CP,164,UX,226,256 

CP,165,UX,227,242,257 

CP,166,UX,228,258 

CP, 1 67,UX,229,243,259 

CP,168,UX,230,260 

CP,169,UX,231,244,261 

CP,170,UX,232,262 

CP,171,UX,233,245,263 

CP,172,UX,234,264 

CP,173,UX,235,246,265 

CP, 1 74,UX,236,266 

CP,175,UX,237,247,267 

CPLIST.ALL 

D,ALL,UY,,,,,ROTX,ROTY ,ROTZ 

D,1,UX 

D,52,UX 

D,73,UX 

D,124,UX 

D,145,UX 

D,196,UX 

D,217,UX 

d,268,ux 

D.289.UX 

D,19,UY 

D,19,UZ 

D,30,UZ 

D,49,UZ 

EP,3 1,5,220 

EP,32,6,220 

EP,33,6,220 

EP,34,6,220 

EP,35,6,220 

EP,3 6,6,220 

EP,37,6,220 

EP,3 8,6,220 

EP,39,6,220 

FINISH 

/solution 

outpr,basic,l 

solve 

finish 
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